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^ \ Let Lq be a closed densely denned symmetric semi-bounded oper- 

i -pH ■ ator with nonzero defect indexes in a separable Hilbert space H. With 

Lq we associate a metric space f2/, that is named a wat>e spectrum 
and constructed from trajectories {u(t)}t>o of a dynamical system 
governed by the equation uu + {Lq)*u = 0. The wave spectrum is 
introduced through a relevant von Neumann operator algebra associ- 
ated with the system. Wave spectra of unitary equivalent operators 
T -j- | are isometric. 

^sO ■ In inverse problems on unknown manifolds, one needs to recover a 

Riemannian manifold Q via dynamical or spectral boundary data. We 
show that for a generic class of manifolds, fi is isometric to the wave 
q ■ spectrum Ql of the minimal Laplacian Lq = — ^\cg°(n\8Ci) acting 

in % = L2(fi), whereas Lq is determined by the inverse data up to 
unitary equivalence. By this, one can recover the manifold by the 
scheme "the data =4> Lq =>• Oi 1& = m $7". 
5_j ■ The wave spectrum is relevant to a wide class of dynamical sys- 

tems, which describe the finite speed wave propagation processes. The 
paper elucidates the operator background of the boundary control 
method (Belishev, 1986) based on relations of inverse problems to 
system and control theory. 
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Introduction 

0.1 Motivation 

The paper introduces the notion of a wave spectrum of a symmetric semi- 
bounded operator in a Hilbert space. The impact comes from inverse prob- 
lems of mathematical physics; the following is one of the motivating ques- 
tions. 

Let Q be a smooth compact Riemannian manifold with the boundary T, 
—A the (scalar) Laplace operator, Lq = — A\c°°(n\r) the minimal Laplacian 
in % = Lz(£i). Assume that we are given with a unitary copy Lq = ULqU* 
in TL = UH (but U is unknown!). To what extent does Lq determine the 
manifold Q? 

So, we have no points, boundaries, tensors, etc, whereas the only thing 
given is an operator Lq in a Hilbert space TL. Provided the operator is 
unitarily equivalent to L , is it possible to extract Q from L ? Such a question 
is an "invariant" version of various setups of dynamical and spectral inverse 
problems on manifolds [2], [I]. 

0.2 Content 

Substantially, the answer is affirmative: for a generic class of manifolds, any 
unitary copy of the minimal Laplacian determines Q up to isometry (Theorem 
1). A wave spectrum is a construction that realizes the determination L =>■ Q 
and, thus, solves inverse problems. In more detail, 

• With a closed densely defined symmetric semi-bounded operator L 
of nonzero defect indexes in a separable Hilbert space % we associate 
a metric space Ql (its wave spectrum). The space consists of the 
so-called eikonal operators (eikonals), so that is a subset of the 
bounded operators algebra whereas the metric on Ql is \\r — 

T> \\<B(H)- 

The eikonals are constructed from the projections on the reachable sets 
of an abstract dynamical system with boundary control governed by the 
evolutionary equation u u + LqU = 0. More precisely, they appear in the 
framework of a von Neumann algebra DT^q associated with the system, 
whereas Ql C 9Tl is a set of the so-called maximal eikonals. The 
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peculiarity is that this algebra is endowed with an additional operation 
that we call a space extension. 

Since the definition of Ql is of invariant character, the spectra VLl 
and of the unitarily equivalent operators L and L turn out to be 
isometric (as metric spaces). So, a wave spectrum is a (hopefully, new) 
unitary invariant of a symmetric semi-bounded operator. 

• A wide generic class of the so-called simple manifolds is introduce^. 
The central Theorem 1 establishes that for a simple Q, the wave spec- 
trum of its minimal Laplacian Lq is isometric to Q. Hence, any unitary 
copy Lq of Lq determines the simple Q up to isometry by the scheme 
Lq =>• =° = ^- I n applications, it is the procedure, which 
recovers manifolds by the BC- method [2] , [4] : the concrete inverse data 
determine a relevant Lq, what enables one to realize the scheme. 

• We discuss one more option: once the wave spectrum of the copy 
L is found, the BC-procedure realizes elements of the space % as 
functions on f2j o 0. Thereafter, one can construct a functional model 
L™ od of the original Laplacian Lq, the model being an operator in 
% raod = L 2 ^(fli o ) related with L through a similarity (gauge trans- 
form). Hopefully, this observation can be driven to a functional model 
of a class of symmetric semi-bounded operators. Presumably, this 
model will be local, i.e., satisfying supp L™ d y C suppy. 



0.3 Comments 

The concept of wave spectrum summarizes rich "experimental material" ac- 
cumulated in inverse problems in the framework of the BC-method, and 
elucidates operator background of the latter. In fact, for the first time Ql 
has appeared in [1] in connection with the M.Kac problem; its later version 
(called a wave model) is presented in [4] (sec. 2.3.4). Owing to its invariant 
nature, Ql promises to be useful for further applications to unsolved inverse 
problems of elasticity theory, electrodynamics, graphs, etc. 

Actually, a wave spectrum is an attribute not of a single operator but 
an algebra with space extension. In the scalar problems on manifolds, this 

Roughly speaking, the simplicity means that the symmetry group of Q is trivial. 
2 In the BC-method, such an option is interpreted as visualization of waves [4]. 
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algebra is commutative, whereas its wave spectrum is identical to Gelfand's 
spectrum of the norm-closed subalgebra generated by eikonals. However, 
it is not clear whether this fact is of general character. The algebras that 
appear in the above mentioned unsolved problems, are noncommutative and 
the relation between their wave and Jacobson's spectra is not understood 
yet. 

By the recent trend in the BC-method, to recover unknown manifolds via 
boundary inverse data is to find spectra of relevant algebras determined by 
the data [5]. We hope for further promotion of this approach. 

1 Wave spectrum 

1.1 Algebra with space extension 

Let Tl be a separable Hilbert space, LatH the lattice of its (closed) subspaces; 
by Pa we denote the (orthogonal) projection onto A G LafH. Also, if A is a 
non-closed lineal set, we put Pa := P c ios^i- By 23(7/) the bounded operator 
algebra is denoted. 

An one-parameter family E = {E l } t >o of the maps E l : Lat"H — > Lat"H is 
said to be a space extension if 

1. E° = id 

2. £'{0} = {0}, t > 

3. t < t' and A C A' imply E l A C E 1 ' A' . 

It is also convenient to regard E l as an operation, which extends the projec- 
tions, and write E 1 Pa = P\ := Pe^a- 

Assume that an extension E is given; let a C LatH be a family of sub- 
spaces. By %l[E, a] we denote the minimal von Neumann operator algebra^, 
which contains all projections {Pa \ A G a} and is closed with respect to E, 
i.e., P G a] implies P t = E l P G *Jl[E,a], t > 0. As is easy to see, 
such an algebra is well defined. As any von Neumann algebra, ^[E, a] is 
determined by the set Proj %l[E, a] of its projections, whereas its additional 
property is 

E l Proj %E, a] C Proj <Jl[E, a], t > 0. 
3 i.e., a unital weekly closed self-adjoint subalgebra of 03(H): see [11] 
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Fix a projection P G %l[E, a]; a positive self-adjoint operator in T-L of the 
form 

t p := / tdP\ 
Jo 

where P* = P*P, is said to be an eikonal^ the set of eikonals is denoted by 
Eik9l[£,a]. 

Let us say that we deal with the bounded case if each eikonal is a bounded 
operator (and, hence, belongs to Ul[E, a}) and the set Eik9T[P, a] is bounded 
in 03(H): 

sup{||r|| | r G EikOI[P,a]} < oo . (1.1) 

Otherwise, the situation is referred to as unbounded case (see the end of sec 
3.4). 

Convention 1 Unless otherwise specified, we deal with the bounded case. 

Recall that the set of self-adjoint operators is partially ordered: for A,Be 
23 ("H) the relation A < B means (Ax,x) < (Bx,x), x G ri. Any monotonic 
bounded sequence A\ < A 2 < . . . , sup ||A, || < oo converges in the strong 
operator topology to s-lim A,- G 23(H) (see, e.g., [6]). 

Extend the name eikonal to all elements of the set s-closEik91[P, a] C 
%l[E, a] and denote the extended set by the same symbol EikOTj-E 1 , a]. An 
eikonal r is said to be maximal if r > r' for any eikonal r' comparable with 
r. Let fl<n[E, a ] C EikOX[£?, a] be the set of maximal eikonals. 

Lemma 1 The set Q<n[E,a] is nonempty. 

Proof By ( II. ip . any totally ordered family of eikonals {r a } has an upper 
bound s-limr a , which is also an eikonal. Hence, the Zorn lemma implies 

The Q3("H)-norm induces the distance \\r — r'\\ in fl<n[E,a] and makes it 
into a metric space, which is determined by the extension E and the initial 
reserve of subspaces a. The space Qyi[E, a ] is the main subject of our paper. 



4 the term is taken from the motivating applications 
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1.2 Extension El 

Here we introduce the space extension associated with a semi-bounded self- 
adjoint operator. Without lack of generality, it is assumed positive definite: 
let 

POO 

L = L*= / A dQx ; (Ly, y) > x \\yf, y G DomL C H , (1.2) 
Jo 

where dQ\ is the spectral measure of L and x is a positive constant. Such 
an operator governs the evolution of a dynamical system 

Va + Lv = h, t > (1.3) 
v\t=o = v t \ t =o = 0, (1.4) 

where h G Ljf ((0, oo); "H) is a "H- valued function of time (control). Its finite 
energy class solution v = v h (t) is represented by the Duhamel formula 

v h {t) = [ L"5 S i n (t- a )L* h(s)ds= (see (TO)) 
Jo L - 1 

* rfs r sinv^-a) ^ ^ ^ £ > Q (L5) 
7o v A 

(see, e.g., [5]). In system theory, v h ( ■ ) is referred to as a trajectory; v h (t) G "H 
is a state at the moment i. In applications, v h describes a wave initiated by 
a source h. Note that in the case of x < the problem (11.31) . (11.41) is also 
well defined but the representation (11.51) is of slightly more complicated form. 
Thus, the assumption x > is accepted just for simplicity. 
Fix a subspace icH; the set 

V\ := {v h {t) | h G 4 oc ((0,oo);^)}, t>0 

of all states produced by ^4-valued controls is called reachable (at the moment 
t). Reachable sets increase as A increases and/or t grows. Indeed, the 
representation (11.51) easily implies 

v T i h (t) = {T $ v h )(t), £>0, (1.6) 

where 7^ is the right shift operator in L!f c ((0, oo); %): 



(%g) (t) 



0, < t < ( 

g(t-£), t>£ 
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with f > 0. For < t < t' and A C ^4', we have 

9 « fc (t) = fa V*) (0 = (see (USD) = v 7 *-**^) G Vj C Vjf, , 
i.e., the inclusion 

V* C VX , < t < t' (1.7) 

does hold. 

Define a family El = {-E*}t>o of the maps E 1 : LafH — >■ LaXH by 

E°A:=A, E l A := closV^, £>0. (1.8) 



Lemma 2 zs a space extension. 



Proof The properties 1 and 2 (see section 1.1) easily follow from the def- 
initions and the obvious relation VL\ = {0}; the property 3 for < t < t' 
is seen from (11.71) . Thus, it remains to verify 3 for t = 0, i.e., to check that 
= t < H and A C A' leads to E°A C E 1 ' A' or, the same, that A C EM 
for all r > and .A 7^ {0}. 

Take a nonzero y £ A and consider (II. 3p . (II. 4p with the control h £ (t) = 
(f e (t) y, where 

t G [0,r-2e) 
t G [r — 2e, r — e) 
t G [r — £, r) 
£ G [r, 00) 



ro, 



v0, 



(e > is small); note that 



<p e {t)f{t)dt -/'(r) 



£->0 



for smooth /'s, i.e., y? £ (£) converges to <5'(t — r) as a distribution. Define 



sin[VA (r - £)] 

7x 



^ e (t) dt 



Z cos (v 7 ^) -cos(v / A2£) - 1 



and note that ip e {X) — >■ 1 uniformly w.r.t. A in any compact segment [x, TV]. 

£— >0 
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Therefore, one has 



\y - v he: (r)\\ = (see(fL5ll) 



dt 



sin[VA (r-t)] 



dQx[ve(t)y] 



ip £ (X)dQ x y 

POO 

/ |i-^(A)| 2 d||g A?/ || 2 -> o 

I e— >0 



[1-V £ (A)] dQ x y 



by the properties of ip £ . The order of integration change is easily justified by 
the Fubini Theorem. 

Thus, y = \imv he (r), whereas v he (r) G E r A holds. By the closeness of 

E r A, we get y%E r A. Hence, A C EM. □ 

So, with each positive definite operator L one associates the certain space 
extension El by (11. 8p . 



1.3 Algebras 9Tl ; x> and 0Tl 



Return to the system ( ll.3p -( }L4|) and fix a nonzero subspace Z> G LaiH that 
we'll call a directional subspace. It determines a class 



A^c := {h G C°° ([0, oo); V) | supp h C (0, oo)} 
of smooth D-valued controls vanishing near s — 0, and the sets 



;i.9) 



U % := yi(t)-v h "(t) 

h(t) — / sin 
Jo 



heM v \ = ( see (JT5) ) 



(t-s)L*] /i"(s)ds 



?ieMp , t>0 (1.10) 



(here (•)' := 4-), which we also call reachable^. As can be easily derived 
from (II. 6p . the sets U l increase as t grows. 



Now, take 



a-L.v '■— {closW } t >o C LatH 



5 The meaning of this definition and term is clarified later on, when we deal with systems 
with boundary control in sec 2.2 
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in capacity of the initial family of subspaces (see sec 1.1). The pair L,T> 
determines the algebra with extension 

m L , v := %E L , a LtT >] 

and its eikonals Eik^t^p. Assuming that fll.ip holds for 9Tl,x>, the set of 
maximal eikonals 

is well defined and called a wave spectrum of the pair L, T>. The operator 

poo 

r 9 := / tdP w (1.11) 



is said to be a boundary eikonal. The set 

dQ LtV := {t E Q L , V \t>t 9 } c n LtV (1.12) 

is a boundary of the wave spectrum. 

Thus, each pair L,V determines an algebra with space extension OT^-p 
and all corresponding attributes. 

Let Lq be a closed densely defined symmetric semi-bounded operator with 
nonzero defect indexes n± = n < oo. As is easy to see, such an operator 
is necessarily unbounded. For the sake of simplicity, it is assumed positive 
definite: (L y,y) > x\\y\\ 2 , y G DomL with x > 0. Let L be the extension 
of L by Friedrichs, so that L = L* > xl and 

L cLcL* (1.13) 

holds [B]. Also, note that 1 < dimKerLg = n < oo. Taking 

V := KerL* 

as a directional subspace, we can constitute the pair L, Ker Lq, which deter- 
mines the algebra 

and its eikonals Eik9lL . If (11.11) holds, the set of maximal eikonals 

^L := &L,Kcr:L* 
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is well defined and referred to as a wave spectrum of the operator Lq; its 
subset 8VLl is a boundary of the wave spectrum. 

So, with every L of the above-mentioned class, one associates the algebra 
9Ti and the algebra eikonals EikOXi . If the latter set is bounded, the 
operator L possesses the wave spectrum Vt Lo ^ 00. 

2 DSBC 

2.1 Green system 

Dynamical systems with boundary control (DSBC) are defined in the next 
sections 2.2 and 2.3; here we introduce a basic ingredient of the definition. 
The ingredient is a collection {H, Q; A, r , Ti} of the separable Hilbert spaces 
% and Q, and the densely defined operators A : H — > H, T k : H — > Q {k = 
0, 1) connected through the Green formula 

(Au, v) H - (u, Av) H = (r w, Tiv) g - (Tiu, T v) g 

(see [IQj). Such a collection is said to be a Green system; Q and T k are 
referred to as a boundary values space (BVS) and the boundary operators 
respectively. In the applications, which we deal with later on, the following 
is also provided: 

1. DomT/j D Domi holds, whereas A is such that the restriction 

^iKerronKerTi =: Lq 

is a densely defined symmetric positive definite operator with nonzero 
defect indexes and A = Lq is valid ('bar' is the operator closure) 

2. the restriction 

coincides with the Friedrichs extension of Lq, so that we have 

L CLCL* Q = A, (2.1) 

whereas L _1 is bounded and defined on % 
6 However, there are examples in applications, in which Q,l consists of a single point. 
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3. for the subspaces A : = Ker A and T> := KerLg, the relations 

clos.4 = £>, closFVl = £ (2.2) 

are valid. 

These properties are in consent with the Green system theory by V.A.Ryzhov, 
which puts them as the basic axioms. Note, that there are a few versions of 
such an axiomatics but the one proposed in [12] is most relevant for applica- 
tions to the forward and inverse multidimensional problems of mathematical 
physics. 

Convention 2 In what follows we deal with the Green systems, which satisfy 
the conditions 1-3. 

As is shown in [12J, the axioms provide the following: 

• the map II := (TiL^ 1 )* : Q — )■ % is bounded, whereas Ran II is dense 
in V 

• the subspace A admits the characterization 

A = {y e Dom A | IW y = y} (2.3) 

• since L is the extension of L by Friedrichs, we have 

DomLo = L~ l [U QV], L = L\ L -i [HeV] (2.4) 

that easily follows from the definition of such an extension (see [6]). 

Example Let Q be a C°°-smooth compact Riemannian manifold with 
the boundary T, A the (scalar) Laplace operator in H := L 2 (fi), v the 
outward normal on T, Q := L 2 (T). Denote 

A = -A\ HHQ) , r :=(.)|r, r i: =|-(-)|r, 

so that r 01 are the trace operators. The collection {H, Q\ A, T , Tx} is a 
Green system, whereas other operators, which appear in the framework of 
Ryzhov's axiomatics, are the following: 

7 H k { . . . ) are the Sobolev classes; Hg(Q) = {y e H 2 (n) \ y = \Vy\ = onT}. 
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is the minimal Laplacian; 

L = — ^\H 2 (n)nH^(n) 
is the self-adjoint Dirichlet Laplacian; 

L = -A\{ y eH\AyeH} 

is the maximal Laplacian; 

A = {yEH 2 (n)\ Ay = 0} 

is the set of harmonic functions of the class H 2 (Q); 

V = {y G % | Ay = 0} 

is the subspace of all harmonic functions in L 2 (f2); II : Q — > % is the harmonic 
continuation operator (the Diriclet problem solver): 

Hip = u : Au = in Q, u\r = ■ 
2.2 Evolutionary DSBC 

Note in advance that the goal of the sections 2.2, 2.3 is not to obtain results 
but motivate the above introduced objects and notions. That is why the 
presentation ignores certain technical details. 

Definition The Green system determines an evolutionary dynamical 
system with boundary control of the form 

u a + Au = inH, < t < oo (2.5) 

u\ t =o = u t \ t =o = in H (2.6) 

T u = f(t) in G, < t < oo, (2.7) 

where / is a boundary control (^-valued function of time), u = u^(t) is the 
solution (wave). 

Assign / to a class J- + if it belongs to C°° ([0, oo); Q), takes the values in 
r DomA, and vanishes near t — 0, i.e., satisfies supp / C (0, oo). Also, note 
that / e F+ implies II (/( ■ )) G Mr>, where V = KerL^ and Mv is defined 

by m>- 
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Lemma 3 For f e T+, the classical solution v? to 'problem l[2.5\) -(2.7) is 
represented in the form 



u f(t) = h{t) 



L 2 sin 



(t-s)L2 



h"(s) ds 



t > 



(2.8) 



with h :=n (/(•)) e M v - 

Proof Introducing a new unknown w = w'(t) := u'(t) — H {f{t)) and taking 
into account ( 12. 3ft . we easily get the system 



i% + Aw = -II (/«(*)) 
u>|t =0 = w t \ t =o = 
T u; = 



in H, < t < oo 
in U 

in ^, < t < oo . 



With regard ro the definition of the operator L (see the axiom 2), this problem 
can be rewritten in the form 



in H, < t < oo 
in U 



w a + Lw = -h tt 
w\ t =o = w t \ t =o = 

and then solved by the Duhamel formula 



w f (t) = - 



Returning back to u* = + II/, we arrive at (12.81) . □ 
The sets 







/ L 2 sin 


\t-s)L*' 


to 





{u f (t)\ fef+} = (see([2lD) 



h(t) - \ L~2 sin (t - s)Lz 



h"(s) ds 



h = Uf(-), fe^ 



t > o 

(2.9) 



are said to be reachable from boundary. In the mean time, the Green system, 
which governs the DSBC, determines the certain pair L, V, which in turn 
determines the family {U 1 } by (11.101) . Comparing the definitions, we easily 
conclude that the inclusion C U l holds. Moreover, the density properties 
(12.21) enable one to derive 



closWl = closW* 



t > 



(2.10) 
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and it is the relation, which inspires the definition ( ll.lOp and motivates the 
use of the term "reachable set" for U 1 in the general case, where neither the 
boundary value space nor the boundary operators are defined. 

Illustration Consider the Example at the end of sec 2.1. The DSBC 
( I2.5p - (l2.7p associated with the Riemannian manifold is governed by the wave 
equation and is of the form 

u tt -Au = Q in fix (0,oo) (2.11) 

u\ t =o = u t \ t=0 = in Q (2.12) 

w|r = f(t) for < t < 00 (2.13) 

with a boundary control / G T := Li, 00 ((0, 00); L 2 (T)); the solution u = 
w(x,t) describes a wave, which is initiated by boundary sources and propa- 
gates from the boundary into the manifold with the speed 1. For / G J-+ := 
C°° ([0, 00); C°°(r)) provided supp/ C (0, 00), the solution is classical. 
By the finiteness of the wave propagation speed, at a moment t the waves 
fill the near-boundary subdomain 

0*[r] := {x G Q\ dist (x,T) < t} . 

Correspondingly, the reachable sets U 1 ^ increase as t grows and the relation 

U\ C H % , t>0 (2.14) 

holds, where H} := clos{?/ G H \ suppy C fi*[T]}. 

So, if the pair L, T> (or the operator Lq) appears in the framework of a 
Green system, then {W*} introduced by the general definition (11.101) may be 
imagine as the sets of waves produced by boundary controls. The question 
arises, what is the meaning of the corresponding wave spectrum Ql,t> (— 
fi^o)? In a sense, it is the question, which this paper is written for. The 
answer (postponed till section 3) is that Ql is a wave guide body, in which 
such waves propagate. 

Controllability Return to the abstract DSBC (I2.5p - (l2.7p and define its 
certain property. The definition is premised with the following observation. 
Since the class of controls J- + satisfies -^■J r + = J-+, the reachable sets ( 12. 9ft 
satisfy AU\ = U\. Indeed, taking / G J-+ we have 

Au f (t) = ( see (E31) ) = -u f tt {t) = u- f "(t) eW{ (2.15) 
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and, by the same relations, 

u f (t) = Au\t) 

with g = — (Jq) 2 / G F+. Hence, the sets reduce the operator A, so that 
its parts A\ u t are well defined. 

The DSBC f!2 .5D — (I2.7P is said to be controllable from boundary for the 
time t = T if the (operator) closure of the part A\ u t coincides with A, i.e., 
the relation 

clos { {u f (T) , Au f (T)} | / G JV} = clos graph A (2.16) 

(the closure in H x W) is valid, where 

graph A := {{y, Ay} \ y G DomA} . 

Controllability means two things. First, since A is densely defined in Tt, 
(12. 161) implies 

closW^ = H, t>T, 

i.e., for large times the reachable sets become rich enough: dense in %. 
Second, the "wave part" A\ u t of the operator A, which governs the evolution 
of the system, represents the operator in substantial: it coincides with A up 
to closure. 

In applications to problems in bounded domains, such a property "ever 
holds" (for large enough times T). In particular, the system (I2.lip - fl2.13p is 
controllable for any T > maxdist (x,T) (see [2], jl]). 

Let us represent the property (1 2 . 1 6 [) in the form available for what follows. 
At first, with regard to (12.11) and (12.151) . it can be written as 



feF+\ = graph L* . (2.17) 



clos l{u f {T),u- f "{T)} 

Further, for each t > 0, introduce a control operator ("input — > state" map) 
W t : J' T-L, Dom W l = J r +, 

W l f := u f (t) . 
In terms of this map, (12.171) takes the form 



clos<i{^ r /,W^(-f)} 



/G^+ = graph L*. (2.18) 



16 



The control operator can be regarded as an operator from a Hilbert space 
J= l := L 2 ((O,t);0) to the space H with Domlf' = J=\ := {f\ m \ f E 
As such, it can be represented in the form of the polar decomposition 

w t = u l \w l \ , 

where |W*| := ((W^W 1 ) 1 and U l is an isometry from closRan|iy*| C J 7 * 
onto closRaniy* C H (see, e.g., jH])- For t = T, one has closRaniy T = H, 
so that U T is a unitary operator from the (sub)space H := closRan \W T \ 
onto "H. Correspondingly, the operator 

LI := (U T YL* Q U T 

acting in H, turns out to be unitarily equivalent to Lq. As result, f)2.18p can 
be written in the final form 



clos \{\W T \f, \W T \(-f")} 



feF + \ = graph L* . (2.19) 



As a consequence, we get 

Proposition 1 // the DSBC \2.!fy -$7fy is controllable for the time T then 
the operator \ W T \ determines the operator Lq up to unitary equivalence. 

Response operator In the DSBC ( I2.5p -( 12T71) . an "input— ^-output" cor- 
respondence is described by the response operator R : T — ^ J-, Domi? = J P + , 

{Rf)(t) := rx(^(t)) , t>0. 

Also, the reduced operators 

R f '■= W) \ [o,t] 

are in use and play the role of the data in dynamical inverse problems. 

The key fact of the BC-method is that the operator R 2t determines the 
operator (W t )*W t through a simple and explicit relation: see [2], [3], [I]. 
Hence, R 2t determines the modulus \W t \. Combining this fact with the 
Proposition (TJ we arrive at 

Proposition 2 // the DSBC $U$-W]§ is controllable for the time T then 



its response operator R 2T determines the operator Lq ( and, hence, the oper- 
ator Lq = Lq* and its Friedrichs extension L) up to unitary equivalence. 
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As illustration, the response operator of the DSBC (12.1ip - (l2.13j) is 



dv 



rx[0,2T] 



By the aforesaid, given for a fixed T > maxdist (x,T) this operator deter- 
^eQ 

mines the operator L up to a unitary equivalence. 
2.3 Stationary DSBC 

Our presentation follows the paper [12J. The basic object is the Green system 
{H, Q; A, T , ri} and the associated operators L , L (see sec 2.1). 

Definition Along with the evolutionary DSBC, one associates with the 
Green system the problem 

(A - zl) u = in H, z E C (2.20) 

T u = (f in £ (2.21) 

that is referred to as a stationary DSBC. For ip G r Dom A and z E C\ spec L, 
such a problem has a unique solution u = v%, which is a Dom A -valued func- 
tion of z. 

Weyl function The "input — > output" correspondence in the system 
fl2.20p - fl2.2ip is realized by an operator-valued function 

M(z)ip := Tiu^ , z specL 

that is called Weyl function and plays the role of the data in frequency 
domain inverse problems. 

The following fact proven in [12] is of crucial value. Recall that a sym- 
metric operator in H is said to be completely non-selfadjoint if there is no 
subspace in H, in which the operator induces a self-adjoint part. 



Proposition 3 // the Green system, which determines the DSBC fiK 



112.21}) , is such that the operator L is completely non-selfadjoint, then the 
Weyl function of the DSBC determines the operator L up to unitary equiv- 
alence. 
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Illustration Consider the Example at the end of sec 2.1. The DSBC 
(12.20p -( l2.2ip associated with the Riemannian manifold is 

(A + z)u = mQ (2.22) 
u| r = <p, (2.23) 

where A = —A\j^(ii)- The operator L = — A\ H 2^ is completely non- 
selfadjoint. Indeed, otherwise there exists a subspace /C C H such that the 
operator Lq := —A\ KnH 2^ ^ O is self-adjoint in fC. In the mean time, 
Lq is a part of L, the latter being a self-adjoint operator with the discrete 
spectrum. Hence, specLg is also purely discrete; each of its eigenf unctions 
satisfies — A0 = A0 in Q and belongs to Hq(Q). The latter implies 0=^ = 
on T, which leads to = by the well-known E.Landis uniqueness theorem 
for solutions to the Cauchy problem for elliptic equations. Hence, Lq = O 
in contradiction to the assumptions. 
The Weyl function of the system is 



M(z)(p 



dv 



(z ^ spec L) 



By the aforesaid, the function M determines the operator L up to a unitary 
equivalence. 

Besides the Weyl function, there is one more kind of inverse data associ- 
ated with to the DSBC (IZ22jl - (IZ251> . Let {X k }^ =1 : < Ai < A 2 < A 3 < 
• • ■ — > oo be the spectrum of the Dirichlet Laplacian L, {0fc}fcli : Lcfy^ = 
Xk4>k its eigenbasis in H normalized by (fa, 4>i) = 5m- The set of pairs 



£n := < A 



dv 



k=l 



is called the (Dirichlet) spectral data of the manifold Q. The well-known fact 
is that these data determine the Weyl function and vice versa (see, e.g., [T2"]). 
Hence, determines the minimal Laplacian L up to unitary equivalence. 
However, such a detrmination can be realized not through M but in more 
explicit way. 

Namely, let U : U -»• U := h, 



Uy = y:={(y,fa)}? =1 
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be the Fourier transform that diagonalizes L: 

L := ULU* = diag{Ai, A 2) ...}. (2.24) 
For any harmonic function a £ A, its Fourier coefficients are 

[a,<p k ) = -— I a——dT 

that can be verified by integration by parts. With regard to the latter, 
the spectral data determine the image A := UA C H and its closure 
T> = UT> = cloSv4, so that the determination 

occurs. In the mean time, ( 12 .4[) implies 

L = U*L U = Zlz-x^ (2.25) 

by isometry of U. Thus, L is a unitary copy of L constructed via the 
spectral data. 

3 Applications 

3.1 Inverse problems 

In inverse problems (IP) for DSBC associated with manifolds, one needs to 
recover the manifold via its boundary inverse data B Namely, 

IP 1: given for a fixed T > maxdist (x, V) the response operator R 2T of the 
system (12. 1 1 [) — ( ]2 . 13j) . to recover the manifold Q 

IP 2: given the Weyl function M of the system f l2.22p -( 12.23p . to recover the 
manifold Q 

IP 3: given the spectral data Eq, to recover the manifold Q. 

The problems are called time-domain, frequency-domain, and spectral respec- 
tively. 

8 In concrete applications (acoustics, geophysics, electrodynamics, etc), these data for- 
malize the measurements implemented at the boundary. 
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Setting the goal to determine an unknown manifold from its boundary 
inverse data, we have to keep in mind the evident nonuiqueness of such a 
determination: all isometric manifolds with the mutual boundary have the 
same data. Therefore, the only reasonable understanding of "to recover" is 
to construct a manifold, which possesses the prescribed data [I] . 

As we saw, the common feature of the problems IP 1-3 is that their data 
determine the minimal Laplacian L Q up to unitary equivalence. By this, each 
kind of data determines the wave spectrum Ql up to isometry. As will be 
shown, for a wide class of manifolds the relation Ql 1S ^? 1 f2 holds. Hence, 
for such manifolds, for solving the IPs it suffices to extract a unitary copy 
Lq from the data, find its wave spectrum f2j '= ^l ; an d thus to get an 
isometric copy of Q. It is the program for the rest of the paper. 

3.2 Simple manifolds 

Recall that we deal with a compact smooth Riemannian manifold Q with the 
boundary T; vol is the volume in Q. Also, recall some definitions. 
For a subset A C Q, denote by 

W[A] := {x G tt \ dist(x,A) < r} 

the metric neighborhood of A of radius r > and put fi°[A] := A. Note 
that whatever A be, its neighborhood is an open set with the zero volume 
boundary: 

vol dtt r [A] =0, r > (3.1) 

[8]. By A^ we denote the set of its interior points: x G A b if there is an e > 
such that fi £ [{x}] C A. For a system a C 2 n , we define a b := {# | A G a}. 

Let Y be a set, H C 2 Y a system of subsets. The system S is said to be 
an algebra if 

• Y G S 

• A,BeE implies Y\A, Af]B G 2 (and hence 0, A U B G E). 

For a family a C 2 y , by S[a] we denote the algebra generated by this family, 
i.e., the minimal algebra that contains a. As is known, H[a] consists of the 
sets of the form U^ =1 D^ =1 A nm , where A mn or F\A mn belong to a (see, e.g., 

0)- 



21 



Return to the manifold. The following is an universal process that asso- 
ciates with Q a certain system ar C 2^ and that we refer to as a Procedure 
1. The process is a consequent repetition of the same operation a that acts 
as follows: for a given family a G 2 Q , 

1. constitute the algebra E[a] and go to the system S b [a] 

2. construct the system of the neighborhoods 



that is the product of the operation a. 

Such an operation is of the following important feature. We say a set A C Q 
to be regular (and write A G 7Z C 2 Q ) if vol A > and vol cM = holds. 
Note that 7Z is an algebra. The feature is that by the definition of the 
minimal algebra and property (I3.ip . the inclusion a <Z7Z implies a[a] C 7Z. 
Also, note that the passage E[a] — > S b [a] removes the zero volume sets, which 
can appear in the algebra 5[a]. 
Now, we describe 

Procedure 1: 

Step 1 Take the family of boundary neighborhoods 7 := {f2*[T]} t >o C 2 Q 
and construct the system a[~f] 
Step 2 Construct a 2 [7] := o"[o"[7]] 
Step 3 Construct cx 3 [7] := er[cr[cr[7]]] 

Final Step Constitute the system 



that is the end product of the Procedure 1. As is easy to see, the constructed 
system consists of regular sets, and is determined by the metric in Q and the 
"shape" of its boundary T. 

A system a C 2 is said to be a net if for any point i6ll there exists a 
sequence {u)j} c *L l C a such that volwj > 0, cu\ D u 2 D . . . , diamco,- — > 0, 
and x G f] Uj. 



We say the manifold f2 to be simple, if the system ar is a net. The 
following is some comments on this definition. 



a[a] := {SI* [A] | A G S b [a], t > 0} 



00 
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The evident obstacle for a manifold to be simple is its symmetries^. For 
a ball Q = {x £ R n | \x\ < 1}, the system ar consists of the sets r) x S 11 ^ 1 , 
where rj C [0, 1] is a sum of positive measure segments. Surely, such a system 
is not a net in the ball. A plane triangle is simple iff its legs are pair-wise 
nonequal. Sufficient and easily checkable conditions on the shape of Q C M. n , 
which provide the simplicity, are proposed in [1J. They are also available for 
Riemannian manifolds and show that simplicity is a generic property: it can 
be reached by arbitrarily small smooth variations of the boundary T. 

For an A C £1, define a distant function 

d>A( x ) '■= dist (x, A) , £ £ . 

By the well-known properties of the distance on a metric space, distant func- 
tions are continuous: (1a £ C(Q). The latter space is a Banach algebra (with 
the sup- norm). Let CV(fi) be the (closed) subalgebra in C(Q) generated by 
the family {d,A \ A £ ar}. The following is the result, which in fact inspires 
the notion of simplicity. 

Lemma 4 // the manifold Q is simple then the equality 

C r (0) = C(ft) (3.2) 

holds. 

Proof For x, x' £ O, i / x', choose uj £ ar such that x £ u and x' £" cJ, 
what is possible since the system ar is a net. One has = cL(x) and 
duj(x') > 0, so that Cr(^) distinguishes points of Q. Therefore, by the Stone 
theorem, the sum CV(fi) V {constants} coincides with C(Q). In the mean 
time, assuming {constants} (jL Cr(^), we have C(Q) = Cr(fi) + {constants} 
and conclude that Cr(Q) is a maximal ideal in C(Q). By the latter, there 
exists a (unique) point x £ Q such that all functions of CV(fi) vanish at x 
(see, e.g., [H]). Evidently, it is not the case since ar is a net. Hence, we 
arrive at (Q. □ 

3.3 Solving IPs 

Here we prove the basic 



9 Presumably, any compact manifold with trivial symmetry group is simple, but it is a 
conjecture. Note that for noncompact manifolds this is not true. 
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Theorem 1 Let Q be a simple manifold, Lq = —A\ H 2^ the minimal Lapla- 
cian, Ql its wave spectrum. There exists an isometry (of metric spaces) i 
that maps Ql onto Q, the relation i(dQL ) = T being valid. 

Proof consists of the parts I — III. 
I. System ar 

We say a system of subspaces E C Lat ft to be an algebra if 

• ft G E 

• A,B eE implies ft Q A, A n B G E (and hence {0}, A V B G S). 

For a family a C Lat ft, by E[a] we denote the algebra generated by this 
family, i.e., the minimal algebra that contains a. As is known, E[a] consists 
of the subspaces of the form U^ =1 V^f =1 A nm , where A mn or H Q A mn belong 
to a. 

Below we present an universal process that associates with Q a certain 
system of subspaces ar C LatH and that we refer to as a Procedure 1. 
The process is a consequent repetition of the same operation a that acts as 
follows: for a given family a C Lat T~L, 

1. constitute the algebra E[a] 

2. construct the family 

a[a) := {E r A \ A G E[a], r > 0} C Lat ft, 

where {i? r } r >o = -El is the space extension determined by the Dirichlet 
Laplacian L D Lo- 

The latter family is the product of the operation a. Such an operation is of 
the following important feature. 

A subspace A C H is called regular (we write A G 7Z C Lat ft) if 

.A = {y G ft | supp y CA} =: ftA 

with an A G so, a regular subspace consists of functions supported on a 
regular set. Now, we invoke a fundamental property of the DSBC (12. 11ft - 
(I2.13P known as a local controllability^, by which 

10 see [4], sec 2.2.3, eqn (2.21). This property is based upon the fundamental Holmgren- 
John- Tataru theorem on uniqueness of continuation of solutions to the wave equation 
(|2.1ip through a noncharacteristic surface: see [2] for detail. 
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• the embedding (12.14ft is dense, i.e., in our current notation, we have 
closW+ = "Hfi*[r], which shows that all closW+ are regular subspaces. 
As a consequence, by (12 . 1 Q[) we conclude that the subspaces closW* are 
also regular: 

closW* = nn^T] CTZ, t>0 (3.3) 

• for any regular subspace HA, the relation 

E t HA = Wl*[A] , t > (3.4) 

holds. 

The above- announced feature of the operation a is the following. By the 
definition of the minimal algebra and property ( 13 ,4p . if a is applied to a 
family a of regular subspaces, then the result a[a] also consists of regular 
subspaces, i.e., a C 1Z implies a[&] C 1Z. 
Now, we describe 

Procedure 1: 

Step 1 Take the family of subspaces 7 := {closW*} 4 > C 1Z corresponding 
to the boundary neighbourhoods fi*[r] (see (I3.3P ) and construct the family 
«r[ 7 ] 

Step 2 Construct a 2 [7] := o"[<r[7]] 
Step 3 Construct <J 3 [t] := o"[<7[<3"[7]]] 

Final Step Constitute the family 

00 

a r := |J(3- J [7] 

that is the end product of the Procedure 1. As is easy to see, the constructed 
family consists of regular subspaces: dp C 1Z. 

The evident duality of the Procedures 1 and 1, which is intentionally 
emphasized by the notation, easily leads to the bijection 

TZ D a r 3 A < — ► HA G d r C K . (3.5) 



II. Eikonals Recall that for a linear set A C H, by P4 we denote the 
projection in % onto clos^4. 
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By the bijection ( 13. 5ft . each projection P4 with ^ G dr is of the form 
P4 = P%a with A G ar- Hence, as an operator in % = /^(fi), P4 multiplies 
functions by the indicator Xa{~) of the set A (cuts off functions on A), what 
we write as P4 = xa- By (13.41) . the projection P l A = E t P A cuts off functions 
on f2*[A], i.e., P\ = Xn*[A]- As result, the operator 

poo 1-00 

r P ^ = / tdP\ = / *dxn*[A] 
Jo Jo 

multiplies functions by the distant function: 

(rp A y) (x) = d A {x) y(x) , x G Q . (3.6) 

Let £oo(fi) be the algebra of Loo-multipliers, which is a von Neumann sub- 
algebra of 03(H). By Cr(^) an d we denote the subalgebras of Cr(fi)- 
and C(Q)- multipliers respectively, both of them being closed w.r.t. the oper- 
ator norm; so, we have £r(^) C C £oo(^) C 53(H). These subalgebras 
are commutative; also, as is well known, is weakly dense in E^ifl) pTTj . 
Recall that 91l is defined in sec 1.3, whereas Lq which we are dealing with, 
is the minimal Laplacian on Q. 

Lemma 5 If the manifold Q is simple then one has 9Tl = £oo(Q). 

Proof (sketch) By construction of ar, all projections P4 with A G dr 
belong to the algebra 9Tl . Whence, the eikonals of the form ( 13.61) belong to 
OTlo, so that the embedding Cr(^) C 91l occurs. 

By simplicity of fi, one has <£r(^) — £(0), what is just a form of writing 
the assertion of LemmalU So, we have C 0T^ . The latter, by the above 
mentioned w-closeness of £(f2) in £,^(^1), implies £00 (^) C 9Tl . 

All projections belonging to £oo(f2) are of the form P = \a with a Borel A 
[TT] . Using the relevant generalization of the relation (13.41) on Borel positive 
volume sets, one can show that the algebra £oo(f2) is invariant w.r.t. the 
extension El- Since 9Tl is a minimal .E^-invariant algebra, the embedding 
£oo(fi) C m Lo yields £,^(£1) = <tt Lo . □ 

Corollary 1 The set of eikonals Eik9T^ consists of the operators of the 
form Ii3. 6\) with the Borel A 's. 

III. Wave spectrum We omit the proof of the following simple fact: 
the eikonal Tp A = d& is maximal iff A is a single point set, i.e., A = {xq} for 
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a Xo G fi. Such an eikonal is denoted by r xo , so that we have 

^L = { r x I g f2} 
by the definition of a wave spectrum, whereas a map 

« : Ql 3 ^ H> x G 

is a bijection. Recall that the distance between eikonals t',t" G Ql is 
||t' — r"|| and show that i is an isometry. 

Take x',x" G f2; the corresponding eikonals t x i,t x " G f2 Lo multiply func- 
tions by g^' and d x » respectively For any x G Q, by the triangular inequality 
one has (^(x) — d x n(x)\ < dist (x',x"); hence for y G "H = L 2 (fi) one easily 
has || [t x i — t x "] y\\ < dist (x', x") ||y|| that implies 

1 1 TV — t x "|| < dist (x',x") . (3.7) 

Choose Uj C f2 such that lo\ D uj 2 D . . . , diamwj — >■ 0, and x' = ^j>i<^j', 
then put y 3 - = llxjll" 1 ^') where Xj is the indicator of Uj. As is easy to see, 
ll^lr'Z/j II ~~ * an d ll r x"l/ill — > dist (x',x") holds as j — > oo that implies 

II [tx> - t x »\ yj\\ -> dist (x',x") , ||%-|| = 1. (3.8) 

Comparing (I3.7P with (13. 8p . we arrive at ||7v — 7V/|| = dist (x',x") and con- 
clude that the map r xo >— >■ xo is an isometry. 

Recall that the boundary of the wave spectrum is introduced in sec 1.3 
by fll.lip and fll . 12j) . In our case, in accordance with (13 .3p the projections 
Put cut off functions on the near-boundary subdomains fi*[T] and, corre- 
spondingly, the boundary eikonal t 9 = J^°tdPut multiplies functions by 
dr := dist(-,T). As is evident, the the relation d XQ > dr holds in Q iff 
Xo G T. Equivalently, the eikonal r Xo G Ql satisfies t Xq > t 9 iff Xo G T. This 
means that the isometry i maps the boundary dQi of the wave spectrum 
onto the boundary T of the manifold. 

Theorem 1 is proved. □ 

Regarding non-simple manifolds, note the following. If the symmetry 
group of Q is nontrivial then, presumably, Ql is isometric to the properly 
metricized set of the group orbits. Such a conjecture is supported by the 
following easily verifiable examples: 
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• for a ball fl = {i 6 K" \x\ < r}, the spectrum Ql is isometric to 
the segment [0,r] C M, whereas its boundary 8Ql is identical to the 
endpoint {0} 

• for an ellipses f2 = {(x,y) G K 2 ^ + w < 1}, fi io is isometric to its 

quarter Q(l{(x,y) \ x > 0, y > 0}, whereas 9f2 Lo ls ^ n | fr + fr = 

1, x > 0, y > 0} 

• let u be a compact plane domain with the smooth boundary, which 
lies in the (open) upper half-plane let Q be a torus in M 3 , which 

appears as result of rotation of to around the x-axis. Then Vti ^ u 

-. isom rs. 

and o\Il = ou. 

Actually, possible lack of simplicity is not an obstacle for solving the 
problems IP 1-3 because their data determine not only a copy of L but sub- 
stantial additional information relevant to the reconstruction of Q. Roughly 
speaking, the matter is as follows. When we deal with each of these prob- 
lems, the boundary T is given. By this, instead of the algebra %1[El, a] = 91^ 
generated by the family a = {closW*}t>o of the sets reachable from the whole 
T (see f ll.lOp and (12. 8ft ) we can invoke the wider algebra %1[El, a'} 3 9t[-Ex, a] 
generated by the much richer family a' = {closW*}t>n iCrC r 3 a of the sets 
reachable from any patch a C T of positive measure lj. As result, although 
the equality 9?l = £oo(^) may be violated by symmetries, the equality 
yt[EL, a'] = £oo(^) always holds, whereas the wave spectrum fly\[E L ,a'] turns 
out to be isometric to Q. The latter is the key fact, which makes the recon- 
struction possible: see [5] for detail. 



3.4 Comments and remarks 

A look at isospectrality As at the end of sec 2.3, let specL = {Afc}^ =1 be 
the spectrum of the Dirichlet Laplacian on Q. The question: "Does specL 
determine Q up to isometry?" is a version of the classical M.Kac's drum 
problem The negative answer is well known (see, e.g., [7]) but, as far 
as we know, the satisfactory description of the set of isospectral manifolds is 
not obtained yet. The following is some observations in concern with such a 
description. 

11 U] J consists of the solutions (waves) u^(t) produced by the boundary controls / sup- 
ported on a x [0, oo) 
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Assume that we deal with a simple Q. In accordance with Theorem 1, 
such a manifold is determined by any unitary copy Lq of the operator Lq C L. 
If the spectrum of L is given, to get such a copy it suffices to possess the 
Fourier image T> = UT> of the harmonic subspace in H = l 2 - see (" 12 . 25[) EE In 
the mean time, as is evident, if Q and Q' are isometric, then the corresponding 
images are identical: T> = T>' . Therefore, Q and Q' are isospectral but not 
isometric iff T> ^ V . In other words, the subspace D is a relevant "index" 
that differs the isospectral manifolds. 

To be a candidate on the role of the harmonic functions image, which is 
admissible for the given L = diag {Ai, A2, • • • } (see (I2.24p ). a subspace V C h. 
has to possess the following properties: 



1. a lineal set : 



h. ev 



is dense in l 2 (see (I2.25P ). whereas 

replacement of T> by any wider subspace T>' D T> leads to the lack of 
density: clos£.g, 7^ / 2 

2. extending an operator L\c- by Friedrichs, one gets L (see ( I2.25P ). 

In the mean time, taking any subspace V C 1% provided 1 and 2 0, one can 
construct a symmetric operator L Q by ( I2.25P and then find its wave spectrum 
as a candidate to be a drum. However, the open question is whether 
such a "drum" is human (is a manifold). 

Model Once again, let Q be simple. By Lemma |5l the algebra 0^ is cyclic, 
i.e., possesses the elements g eH = L 2 (fl) such that 

clos {Pg \ Pe ProjtftiJ = U. 

This enables one to realize elements of H as functions on the wave spectrum 
of 9Tl by the following scheme: 

• Fix a cyclic g G % and endow Ql with a measure fi as follows. For 
a maximal eikonal r = J^°tdP* £ Ql and a ball S r [r] := {r' £ 
^L I Ik' - < r} put 

/ i( J B r [r]):=(P T ^,^ (3.9) 



12 It is the fact, which is exploited in [I] 

13 such subspaces do exist (M.M.Faddeev, private communication) 
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and then extend /i to the Borel subsets of Ql . As is easy to check, the 
equality 

fji(B r [r]) = / \g\ 2 dvol (3.10) 

holds, where r and x are related through the isometry i : Ql — > Q 
established by Theorem [TJ x = 

So, we have a model space "H mo d := L 2 ^(flL )- 
For a y E H, define its image Iy E 7i mo d by 

(Jy)(r):=lim^Mk. (3 .H) 

The relations f 13 . 9 j) and (I3.10P easily imply 

f 0rr , u q dvol 

(Iy) (r) = hm . = (iT 1 *) (*) (3.12) 

and, hence, the image map / is a unitary operator from % onto % m od- 



An operator 



L™ od := /Lo/* 



can be regarded as a functional model of the operator L on its wave spectrum 
that we call a wave model. In the case under consideration, we have L™ od = 
g" x L g, i.e., this model is just a gauge transform of the original. As such, it 
is a local operator: 

supp L™ A w C supp w (3.13) 

holds for w E DomL™ od . 
Conjectures 

1. We suggest and hope that "H mo d and L™ od do exist for a wide class 
of symmetric semi-bounded operators L , the locality property (13.131) being 
held. In contrast to the known models (see, e.g., p3]), this one has good 
chances to be of the real use for applications. 

The principal point (and difficulty) is to attach the invariant meaning 
to the limit (13. lip . Presumably, it can be done in the framework of the 
representation 

m Lo = © / 3MMr) ( 3 - 14 ) 
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in the form of a 1-st kind von Neumann algebra. Such a representation diag- 
onalizes 9Tl and is expected to be valid for Lo's coming from mathematical 
physics. As an encouraging example, the Maxwell system in electrodynamics 
can be mentioned: see [I]. Note that Maxwell's 9Tl is noncommutative. 

2. A question of independent interest is whether any von Neumann al- 
gebra with space extension (in particular, DTl ) is of the form (13.14)) . Also, 
in addition to the metric ||r — t'\\, it is reasonable to look for more subtle 
structures on Ql like tangent spaces, differentiate structure, etc. 

3. One more attractive option is to construct a wave model of the abstract 
Green system satisfying Ryzhov's axiomatics, by the scheme 

{%, g ; a, t , rj l n Lo => {n mod , g-, A mod , r- od , rr d } 

with "Hmod = © fo H T duir) and a local A mod . An intriguing point is that 

the boundary 8D,l q is well defined, so that there is a chance to realize T™° d 
as the "true" trace operators. Such a model would provide a canonical real- 
ization of the original system, the realization being determined by its Weyl 
function M and, hence, relevant to inverse problems. 

Unbounded case Return to the sec 1.1. If the assumption (II. ip is cancelled, 
the set Eik^f-E, a] is still well defined but unbounded. In particular, the 
case that all eikonals r = J Q °° t dP t are unbounded operators is realized in 
applications: for instance, it holds if we deal with a noncompact simple 
manifold Q. As result, even though in the mentioned example the relevant 
maximal eikonals do exist, in the general situation we have to correct the 
definition of the wave spectrum £l<n[E, a ]- A possible way out is to deal with 
the regularized eikonals r = J °° j^— dP l with a fixed a > and thus reduce 
the situation to the bounded case. 

A bit of philosophy In applications, the external observer pursues the goal 
to recover a manifold Q via measurements at the boundary T. The observer 
prospects Q with waves u* produced by boundary controls. These waves 
propagate into the manifold, interact with its inner structure and accumulate 
information about the latter. The result of interaction is also recorded at T . 
The observer has to extract the information from the recorded. 

By the rule of game in IPs, the manifold itself is unreachable in principle. 
Therefore, the only thing the observer can hope for, is to construct from the 
measurements an image of fl possibly resembling the original. By the same 
rule, the only admissible material for constructing is the waves w . 
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To be properly formalized, such a look at the problem needs two things: 

• an object that codes exhausting information about Q and, in the mean 
time, is determined by the measurements 

• a mechanism that decodes this information. 

Resuming our paper, the first is the minimal Laplacian L , whereas to decode 
information is to determine its wave spectrum constructed from the waves 
u . It is £Il , which is a relevant image of ft. 

The current paper develops an algebraic trend in the BC-method [5], by 
which to solve IPs is to find spectra of relevant algebras. An attempt to 
use this philosophy for solving new problems would be quite reasonable. An 
encouraging fact is that in all above-mentioned unsolved IPs of anisotropic 
elasticity and electrodynamics, graphs with cycles, etc, the wave spectrum 
Ql does exist. However, to recognize how it looks like and verify (if true!) 

that Ql =° f2 is difficult in view of very complicated structure of the cor- 
responding reachable sets U l . 
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